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Abstract: Quark-gluon plasma during its initial phase after its production in heavy- ion 
collisions is expected to have substantial pressure anisotropies. In order to model this 
situation by a strongly coupled M = 4 super- Yang-Mills plasma with fixed anisotropy by 
means of AdS/CFT duality, two models have been discussed in the literature. Janik and 
Witaszczyk have considered a geometry involving a comparatively benign naked singularity, 
while more recently Mateos and Trancanelli have used a regular geometry involving a 
nontrivial axion field dual to a parity-odd deformation of the gauge theory by a spatially 
varying 9 parameter. We study the (rather different) implications of these two models on 
the heavy-quark potential as well as jet quenching and compare their respective predictions 
with those of weakly coupled anisotropic plasmas. 
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1. Introduction 

In ultrarelativistic heavy-ion collisions, quark-gluon plasma is produced far from equilib- 
rium with strong anisotropy caused by the fact that initially the system expands mainly 
along the collision axis. This complicates enormously any theoretical analysis and makes 
it difficult to decide whether the strong collectivity observed is indeed proving that the 
quark-gluon plasma behaves as a intrinsically strongly coupled, near-perfect fluid and that 
a description based on a perturbative plasma with collective effects from strong gluon fields 
can be ruled out [1]. 

At weak coupling, an anisotropic quark-gluon plasma leads to (non-Abelian) plasma 
instabilities [2, 3, 4, 5, 6, 7, 8, 9] leading to nonperturbatively large fields and turbulent 
behavior [10, 11, 12] which could be responsible for the strong collectivity. Using the 
framework of hard- loop effective theory [13] , experimental signatures such as anisotropic 
photon and dilepton emission [14], momentum broadening of jets [15, 16] as well as heavy- 
quark potentials [17, 18] and quarkonium dissociation [19] have been studied with fixed 
anisotropy as an approximation to the actual dynamical situation. 

At strong coupling, for which conventional lattice gauge theory is of no help for ana- 
lyzing strongly nonequilibrium dynamics, holographic gauge/gravity duality [20, 21] offers 
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the prospect of providing qualitative and semi-quantitative insights. In this framework, 
heavy-ion collisions and the subsequent anisotropic dynamics and eventual thermalization 
of a (super-)Yang-Mills plasma can be modeled e.g. by collisions of shock waves in anti-de 
Sitter space and horizon formation [22]. 

At least for certain observables, it may be meaningful to consider approximations with 
a temporarily fixed anisotropy in the gravity dual. One such attempt to model the effects of 
a system with anisotropic pressures was proposed by Janik and Witaszczyk in [23], where 
the gravity dual involves a comparatively benign naked singularity. In [24] we have studied 
electromagnetic signatures of this model which are qualitatively similar to weak-coupling 
results at high frequencies. 

More recently, Mateos and Trancanelli [25, 26] have proposed a regular gravity dual 
of an anisotropic but equilibrium N = 4 super- Yang-Mills plasma where a stationary 
anisotropy is introduced through a (parity- violating) deformation of the gauge theory with 
a 9 parameter that depends linearly on one of the spatial coordinates. In contrast to the 
singular gravity dual of Ref. [23], this model has a hydrodynamical limit, and in [27] we 
have shown that it involves the remarkable feature that certain components of the viscosity 
tensor break the usual holographic bound of Einstein gravity duals. This model has been 
explored in [28, 29, 30, 31] with regard to the heavy quark potential, the drag force on 
quarks, and jet quenching. (Other holographic models of anisotropic fluids have been 
introduced and studied in Refs. [32, 33, 34].) 

The aim of the present paper is to compare the effect of anisotropy on the heavy quark 
potential and jet quenching in the two holographic models by Janik and Witaszczyk (JW) 
and Mateos and Trancanelli (MT) based on a singular gravity dual and on axion-dilaton 
gravity, respectively. These results are moreover compared with perturbative results for a 
weakly coupled anisotropic plasma for both prolate and oblate anisotropies. 

The organization of this paper is as follows. After briefly reviewing in Sect. 2 the 
two holographic models for an anisotropic super- Yang-Mills plasma of Refs. [23, 25, 26], 
in Sect. 3 we recall the effects of anisotropies on a weakly coupled plasma as described 
by hard anisotropic loop effective theory and we then consider a zero-coupling version of 
the MT model and its consequences. In Sect. 4 we consider the heavy-quark potential in 
the two holographic models as well as in the two weak-coupling models of an anisotropic 
plasma. In the case of the hard anisotropic loops we reproduce results for the real part of 
the heavy quark potential obtained previously in Refs. [17, 18]. In addition we note that 
in certain directions and at large distances the anisotropic heavy quark potential involves 
oscillatory behavior instead of an exponential tail, which is caused by electric plasma 
instabilities. Those are absent in the ^-deformed weakly coupled gauge theory, which 
however also has a nonmonotonic behavior at large distances in the anisotropy direction. 
The holographic models first of all show a completely different large-distance behavior 
since there is a finite distance where the string connecting the heavy quark breaks. Below 
the string-breaking distance, the heavy quark potential in the JW model agrees with the 
hard anisotropic loop results in that for oblate anisotropies the binding is stronger along the 
anisotropy direction than transverse to it and also in that this feature is reversed for prolate 
anisotropy. However, in the MT model as well as in its zero-coupling version the potential 
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is always deeper in transverse directions, for both prolate and oblate anisotropy. In Sect. 5 
we compare the two holographic models with regard to jet quenching and momentum 
broadening parameters. Again we find that in the JW model the effects of anisotropy are 
opposite for prolate and oblate anisotropies, while the MT model is more uniform in this 
respect. However, in the case of anisotropic jet quenching parameters, the weak coupling 
results of Ref. [15, 16] show the opposite trend than those of the JW model while they 
happen to agree qualitatively for an oblate anisotropy in the case of the MT model. Sect. 6 
contains our conclusions. 

2. Two holographic models of strongly coupled anisotropic plasma 

A primary measure of the anisotropy of the boundary field theory is the pressure anisotropy 

A=^-l, (2.1) 

where A > (A < 0) corresponds to an oblate (prolate) plasma. In the following we shall 
recapitulate the main features of the two gravity duals we want to consider and also discuss 
the behavior of the stress-energy tensor in either case. 

2.1 JW model: singular gravity dual 

The metric of the dual geometry given in Fefferman-Graham coordinates 

dg 2 = 7nu(x CT ,u)dxVdx u + du 2 
u 2 

can be related to the expectation value of the stress-energy tensor. Here u is the holographic 
coordinate with u = defining the AdS boundary. The correspondence is such that near 
u = the metric should be given by 

1^,u) = v, + u 4 7 (iV) + 0(u e ). (2.3) 

with 

AT 2 

In [23] Janik and Witaszczyk fixed the (traceless) stress-energy tensor 

(T^(:r CT )> = diag(e, P ± , P ± , P z ) (2.5) 

and therefore the boundary conditions for the Einstein equations with a negative cosmo- 
logical constant, without adding any further matter fields. The most general form of the 
metric respecting the symmetries of the stress-energy tensor is 

ds 2 = -2 (-a(u)dt 2 + c(u) (dx 2 + dy 2 ) + b(u)dz 2 + du 2 ^j . (2.6) 

Solving for the unknown functions a(u), b(u) and c(u) one finds 

a(u) = (l + ^V) 1 / 2 -^ 3 ^/ 4 ^- ^4)1/2+^36^/4 

b(u) = (1 + ^2^)1/2-^/3+^36^/12^ _ A 2 u 4)l/2+B/3-V36=W/12 ^ 
c {u) = (1 + J 42 u 4)l/2+B/6+V36-2B2/12( 1 _ ^2^4)1/2-5/6-^36-252/12^ 
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where A and B are related to the energy density and the pressures by 

A is a dimensionful parameter which in the isotropic case (B = 0) is related to tem- 
perature T according to A = n 2 T 2 /2. Nonvanishing values of the dimensionless parameter 
B characterize the anisotropy of the system, since P z and P± depend differently on B. For 
negative B the plasma is prolate, while it is oblate for positive B. In the following we shall 
consider in particular the values B = y/2, where P z = (A = oo), and B = where 
P± = (A = —1). In a plasma made of free particles, such values correspond to maximal 
anisotropies, but the above geometry permits also negative values of pressure components 
for larger B (limited only by \B\ < \/l8). 

The gravity dual we just described is pathological in the sense that a naked singularity 
appears whenever B does not vanish. The singularity is however benign in the sense 
that one can still choose infalling boundary conditions at the singularity, such that the 
calculation of for example retarded current-current correlators is still possible [24]. It 
turns out that these current-current correlators show no hydrodynamical behavior and 
strongly deviate from the isotropic result even for arbitrarily small -B's for sufficiently small 
frequencies 1 . However, for larger frequencies the results are well behaved and smoothly 
approach the isotropic result for decreasing \B\. In [24] we therefore assumed that this 
background can approximately describe the non-equilibrium physics on short enough time 
scales. Eventually the plasma will evolve towards equilibrium and therefore stationarity 
is not a valid assumption anymore. Asking questions about zero frequency limits, which 
strictly speaking probe an infinite time span are not meaningful in this model. In principle, 
this is also a problem for the potential between two heavy quarks in the plasma which we 
shall compute below, but we expect that the stationary approximation can still provide 
some qualitative insight. 



(2.8) 
(2.9) 
(2.10) 



2.2 MT model: axion-dilaton-gravity dual 

In [25, 26] Mateos and Trancanelli presented a completely regular and well behaved gravity 
dual to an anisotropic but static plasma. Their model is based on the spatially anisotropic 
duals of Lifshitz-like fixed point of [36], but with AdS boundary conditions. This provides 
an anisotropic version of an N = 4 super- Yang-Mills plasma where the anisotropy is kept 
fixed by a parity violating deformation of the gauge theory 

Sgauge = S N=A + 8S, SS = / 6(z)Tv F A F, (2.11) 



1 A similar behavior has in fact been found recently for the out-of-equilibrium production rate of dileptons 
modeled by collapsing shells in AdS [35] . 
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Table 1: Brane set up. 



t= — + — = x + ie-^. (2.12) 



where 6{z) = 2imr>'jz depends linearly on one spatial dimension. riD7 is a constant with 
dimensions of energy and can be interpreted as a density of D7 branes. The complexified 
coupling constant of the SYM theory is related to the axion-dilaton of type IIB supergravity 
by 

8 Airi 
2n 9ym 

For the deformation in (2.11) the axion field x = az ls position dependent with a = 
Xhdt/^Nc written in terms of the 't Hooft coupling A = g\ M N c . Since the axion is 
magnetically sourced by D7 branes, the solution can be considered as a number of such 
branes dissolved in the geometry. It is important to note that since the D7 branes do not 
extend in the holographic direction, they do not reach the AdS boundary and therefore 
do not introduce new degrees of freedom in the N = 4 SYM theory. The brane setup is 
summarized in Table 1. 

The solution for the ten-dimensional bulk geometry that we are eventually interested in 
is a direct product of a five-dimensional manifold M with a negative cosmological constant 
A = — 6/L 2 and S 5 with radius L given by L 4 = 4irg s N c lg. Therefore it suffices to consider 
the five dimensional axion-dilaton-gravity action 

where k 2 = SitG = 4ir 2 /N 2 and L = 1. The line element in the string frame is of the form 2 

1 / fll 1^ \ 

ds 2 = -2 ( - -F(u)£(u)cft 2 + dx 2 + dy 2 + U(u)dz 2 + ^— ) . (2.14) 

In the following we will stop to write the dependence on the holographic coordinate u 
explicitly. We note that reparametrization invariance is already used to fix the coefficient 
in front of dx 2 and dy 2 and that B cannot be set to unity in general. If % = 1 we would 
get an isotropic solution. T is the blackening factor that must vanish at the position of the 
horizon u = Uh- It turns out that all the functions B, J- and T~L can be written in terms of 
the dilaton <p, which itself has to satisfy a third order nonlinear differential equation in u. 

We note that the temperature and the entropy density are well defined since the 
solution under consideration is static. The temperature can be found from the regularity 
condition on the metric after Euclidean continuation and is given by T = \J : ' / (u^ l )\^/Bh/4^r. 
The entropy density is given by a quarter of the horizon area over spatial volume. 

The thermodynamics of this setup is discussed in detail in [26]. To summarize some 
of the most important points: 

2 Here u does not correspond to the holographic coordinate in Fefferman-Graham coordinates as before. 
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Table 2: Choice of parameters and corresponding thermodynamic quantities for approximately 
constant temperature in the anisotropic axion-dilaton gravity dual. 
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Table 3: Choice of parameters and corresponding thermodynamic quantities for approximately 
constant entropy density in the anisotropic axion-dilaton gravity dual. 

• Holographic renormalization brings in a reference scale \x and therefore the stress- 
energy tensor of the boundary theory shows a conformal anomaly (T{t) oc a 4 . The 
energy density and the pressures depend separately onT/n and a/fi. 

• When we keep the temperature constant and increase the anisotropy parameter from 
the isotropic limit a = the pressure anisotropy first always becomes oblate. The 
maximal value of A depends crucially on the temperature. After this initial oblate 
phase there always exists a special value for a where the pressures in transverse and 
longitudinal direction coincide (without the bulk geometry becoming isotropic) and 
if a is increased further the plasma becomes increasingly prolate. However regardless 
of the pressure anisotropy, the bulk geometry of the MT model is uniformly prolate, 
H>1. 

• For small values of a the plasma is unstable against filamentation along the z- 
direction. It is thermodynamically favorable to have regions in z (but infinitely 
extended in x- and y-direction) that are isotropic and regions with a larger value of 
a. However, the interval of a's for which these filamentation instabilities are present 
is smaller than the interval for the oblate plasma. In other words, the prolate phase 
is always stable, but there also exist oblate and stable phases. 

In the following sections we will compute the heavy quark potential and and the jet 
quenching parameter at constant temperature and at constant entropy density. In Table 2 
and Table 3 we present our choice of parameters for these two situations and some related 
thermodynamic quantities (in units with \i = l) 3 . These parameters are chosen so as to 
be in the same ballpark as those considered in Ref. [25, 26]. With different parameters, 

3 The parameter tfi h is related to the dilaton and the anisotropy parameter by 4>h = 4>( u h) + f log a. 
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Table 4: Anisotropy parameters in the hard anisotropic loop effective theory 

also larger positive values of A are possible, however the following results do not change 
qualitatively. 

3. Two weak coupling models of a stationary anisotropic plasma 
3.1 Hard anisotropic loop effective theory 

In a weakly coupled (nearly collisionless) ultrarelativistic gauge theory plasma there is a 
hierarchy of scales, with hard scales p defined as typical energies and momenta of plasma 
constituents, and soft scales gp, with coupling constant g <C 1, pertaining to leading-order 
collective phenomena such as Debye screening and plasmon masses. In thermal equilibrium, 
the effective theory of soft scales is provided by the "hard thermal loop" effective action 
[37]. With anisotropic distribution functions for hard particles, the corresponding "hard 
anisotropic loop" effective theory [13] involves a rich spectrum of stable and unstable 
modes at momentum scales gp, which have been worked out completely for axisymmetric 
deformations of distribution functions of the form [3] 



f(p) = iV/isotv^ + fri/Phard) 



(3.1) 



with anisotropy direction z and some normalization factor N(£) with iV(0) = 1 for vanishing 
anisotropy parameter £. A prolate momentum distribution is obtained for — 1 < £ < 0, 
whereas £ > parametrizes oblate momentum distributions. 

While the pressure anisotropy A is directly determined by £ (see Table 4) , a comparison 
of quantities at different anisotropy is rather ambiguous [38]. This could be done, e.g., by 
keeping the number density or the energy density fixed, but in both cases it also depends 
on whether this is done by adjusting the normalization N or the parameter phard- Keeping 
number densities of hard particles fixed by adjusting N, as done in Ref. [18], leads to 
jV( n )(£) = y/1 + £, whereas constant energy density in hard particles requires N^ e \^) = 

n- 1 ^) with 



no 



'![(! + O" 1 + r 1/2 arctan(Ve)] for £ > 

0" 1 + (-£)- 1/2 atanh(V^)] for £ < 



(3.2) 



Alternatively, one could compare isotropic and anisotropic plasmas by fixing N = 1 
and adjusting phard- Keeping the number density constant requires 



*& = (1 + 1/6 T, 



(3.3) 
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whereas for constant energy density one has 



p£Ld = ^~ 1/4 (e)T, (3.4) 

with T = Phard|^=o- Following Refs. [19, 38], we shall mainly consider the option of rescaling 
Phard in Sect. 4.2. 

At leading order, an anisotropic distribution function of the form (3.1) gives rise to a 
polarization tensor of the form U llu (k) = m 2 D ti pv (uj /\k\^ \k z \/\k\,£) with four independent 
dimensionless structure functions and m 2 D the isotropic Debye mass. 

For any nonzero £, the (chromo-)magnetostatic propagator turns out to involve space- 
like poles and unstable modes corresponding to a filamentation (or Weibel) instability. Be- 
low we shall be interested also in the (chromo-)electrostatic potential given by the Fourier 
transform of the electrostatic propagator 

k 2 + m 2 + m 2 

D 00 (uj = Q,k) = — % ^— (3.5) 

(k 2 + m 2 ^ + m 2 )(k 2 + m 2 ) — m\ 

where m 2 a ^ s are elementary but rather unwieldy functions of £ and k 2 / k 2 (for explicit 
expressions see [17]). This electrostatic propagator has poles at purely imaginary wave 
vectors, corresponding to Debye screening, but additionally poles at a range of real wave 
vectors, which correspond to electric plasma instabilities. The latter arise for wave vectors 
within (outside) a 45° cone about the direction of anisotropy for the oblate (prolate) case. 

3.2 Anisotropic Chern-Simons deformation of weakly coupled gauge theories 

As an alternative model of a plasma with fixed anisotropy we consider the zero-coupling 
limit of a gauge theory with the deformation (2.11) present in the MT model, i.e. a gauge 
theory with Lagrangian 

C = -\F« u F a ^ - -^A^. (3.6) 

This is similar to axion electrodynamics with constant spacelike axion gradient [39, 40] 4 
but in the spirit of the MT model we interpret a, which has dimension of inverse length, 
thermodynamically as a density of some conserved charge distributed along the spatial z 
direction. 

The Lagrangian (3.6) implies an anisotropic dispersion law for gauge boson modes with 
two gauge-invariant branches [40, 42] 



Since u± > 0, there are no tachyonic modes, in contrast to the case of a timelike axion 
gradient [41] and also in contrast to the hard anisotropic loop effective theory. 



4 With timelike axion gradient, this model is known as the Carroll-Field-Jackiw model [41] of a Lorentz 
and parity violating (but isotropic) modification of electrodynamics. 



- 8 - 



Introducing a finite temperature, the free energy per gauge boson turns out to be given 

by 

in the limit T S> a [42]. In the high-temperature limit we can ignore renormalization 
ambiguities, which appear in the 0(a 4 ) terms in the present case of zero coupling (a more 
complete discussion will be given in [42]). With the interpretation of a as a density of a 
conserved charge (as in [26]) which can be increased by compressing the volume such that 
a scales inversely to its longitudinal extent, we obtain the pressure anisotropy 

pf) _ pf) = a 8T = ^__^I + 0[al) (39) 

which indicates a prolate pressure anisotropy at sufficiently large T/a. This is to be con- 
trasted with the holographic infinite-coupling result of the MT model [26] 

P z -P ± = - a ^ + (a 4 ) (3.10) 

which corresponds to an oblate anisotropy when a <C T. 5 

In the following we shall compare the various models in particular with regard to 
anisotropics in the heavy quark potential. In the weak-coupling model given by (3.6), this 
is obtained from the anisotropic electrostatic propagator which has the simple form 

D 00 (oo = 0,k) = ^ 1 — (3.11) 

K ' P + a^l-fcf/P) V ' 

This corresponds to anisotropic Debye screening, without the complication of electric in- 
stabilities which are present in the electrostatic propagator of the hard anisotropic loop 
effective theory. (In fact, there are also no magnetic instabilities in the a-deformed mag- 
netostatic propagator, although as will be discussed in [42] the phase diagram in the a-T- 
plane has regions of metastability and absolute instability against filamentation towards 
inhomogeneous densities a along the direction of anisotropy.) 



4. Heavy Quark Potential 
4.1 Holographic calculations 

We begin by discussing the heavy quark potential obtained from the Wilson-Polyakov 
loop which is dual to a fundamental string with spacelike separated endpoints at the AdS 
boundary [43, 44, 45] for the generic form of a metric describing stationary but spatially 
anisotropic geometries 

ds 2 = g tt (u) 2 dt 2 + g xx {u) (dx 2 + dy 2 ) + g zz (u)dz 2 + g uu (u)du 2 . (4.1) 

5 Outsidc of the high-temperature limit, oblate as well as prolate phases appear, both in the strong- 
coupling MT model and in the zero-coupling case [42]. 
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The action for the hanging string is 



h l s "^ (4 - 2 > 



S = 

with hob = gABd a X M dbX N being the induced metric on the worldsheet. Here the indices 
a, b are either or 1 and M, N = {t, x, y, z, u}. Due to the symmetry in the transverse plane 
we can always choose a coordinate system such that the y-coordinate vanishes. Parametriz- 
ing the string worldsheet by t and u and making a stationary ansatz for x = x(u) and 
z = z(u), we obtain 

S = -j i - 1 J dt du C(x'(u),z'(u),u) (4.3) 



T 



2ira' 

where primes denote derivatives with respect to the holographic coordinate u and T is a 
constant coming from the time integration. We need to find the string profile and therefore 
evaluate the equations of motion for x(u) and z(u), which are of the form 

-gtt{u)gxx{u)x'(u) =U x jC(x'(u),z'(u),u), (4.4) 
-9tt(u)g zz {u)z'(u) = U z C(x'(u),z'(u),u), (4.5) 

Tl x and H z being constants of motion. Disentangling the above equations we end up with 

= nlguu(u)gtt(u)g zz (u) ^ 

9xx{u) {gtt{u)g xx {u) + n|) (g tt (u)g zz (u) + n^) - n^n^ 

z ,2^ = n 2 z guu(u)g tt (u)g xx (u) ^ ^ 

g zz {u) [(gtt(u)g xx (u) + Ul) (g tt (u)g zz (u) + U*) - U 2 x uf 

For a hanging string that connects two spatially separated points at the boundary we 
expect x' 2 (u) and z' 2 (u) to become negative for u > uq. Since the numerator is manifestly 
positive (note that in our conventions gtt(u) is negative for Lorentzian signature) the de- 
nominator has to vanish at some point uq and then becomes negative for increasing values 
of u. This is also in line with the requirement that du/dx = du/dz = at the turning 
point uq. Evaluating the zero in the common factor of the denominators eventually leads 
to a equation that can be written as 

n 2 n 2 

+ = -ga(u ) > 0. (4.8) 



gxx 

{uq) g zz (u ) 
This is the defining equation of an ellipse and therefore 



n x = -gtt(u )g X x{uo) sin 2 cf>, (4.9) 
n z = -gtt{u )g zz (u ) cos 2 4>. (4.10) 



- 10 - 



x'(u) and z'(u) are completely determined by uq and the angle <f>. To obtain the 
functions x(u) and z(u) we can make one further choice, namely that both x and z vanish 
at the turning point uq. It is then easy to find the distance between the two string endpoints 

L = 2^/x 2 (0) + z 2 (0) (4.11) 

and the energy of the configuration 

E re9 - = -- , / du^-g tt (u)g uu (u). (4.12) 

To calculate the action above we integrate from u = to the turning point u = uq 
in (4.4), which covers only half of the string and therefore we have to multiply by two in 
order to obtain the full result. The last term above is the energy of two straight strings 
hanging from the boundary to the horizon at Uh and is necessary to regularize the amount 
of energy of the hanging string. This also means that the connected configuration is 
energetically favored as long as E reg ' < 0. It can be checked easily that for an isotropic 
geometry with g zz (u) = g X x( u ) we recover the already well known expression for the heavy 
quark static potential. If in the anisotropic case we restrict to the simpler cases where the 
string endpoints are either separated exactly along the z or x direction the above equations 
simplify and we reproduce the same solutions as given previously in [28]. Our expressions 
above however are valid for any separation of the string endpoints in the xz-plane. The 
generic situation allows us to probe the geometry by letting the string hang down in the 
bulk and study how it deforms function of u. 

We start by discussing the results for the JW model, the singular anisotropic gravity 
dual. In Fig. 1 we plot the potential between the two heavy quarks, where we have adjusted 
the parameter A of the model such that the energy density is kept constant for different 
anisotropics. Full (dashed) lines correspond to quarks separated along (transverse to) the 
direction of anisotropy. 

We note that in the oblate phase quarks separated along a transverse direction have a 
slightly shallower potential and consequently a smaller dissociation distance. (By dissocia- 
tion distance we are referring to the maximal distance between two quarks, for which it is 
still energetically favorable to be connected by a hanging string in the bulk. 6 ) For prolate 
plasmas the heavy quark potential is instead shallower for longitudinally separated quarks 
than for transverse separations. Evidently the anisotropy only mildly influences the heavy 
quark potential even though we are considering extreme anisotropic plasmas with P z = 
(B = y/2) and P ± = {B = -y/G). In Fig. 2 we have made these small effects more 
conspicuous by plotting the difference in the separation of two quarks at a given potential 
energy compared to the isotropic case. 

Let us finally study the profile of the hanging string in the singular anisotropic geometry 
of the JW model in more detail. Due to the deformation of the spacetime as we go away 

6 Strictly speaking, the dissociation of heavy quarkonia in a medium also depends on the imaginary part 
of the static potential in the real-time formalism which leads to a finite thermal decay width [46, 19] and 
which we ignore by only studying Wilson loops in the Euclidean time direction. 
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Figure 1: Potential energy of heavy quarks 
in the JW model for plasmas with different 
anisotropics but constant energy density. The 
isotropic case corresponds to B = (black 
lines) , oblate anisotropy with P z = to B = v2 
(blue) and prolate anisotropy with P± = to 
B = — v6 (red). Full (dashed) lines refer to a 
separation of the quarks along (perpendicular 
to) the direction of the anisotropy. 



Figure 2: Difference in the distance between 
two connected quarks at a given potential with 
its isotropic value for a plasma with oblate 
(P z = 0, B = \/2, blue) and prolate (P± = 0, 
B = — V6, red) anisotropy at the same energy 
density. Full (dashed) lines correspond to a sep- 
aration of the quarks along (perpendicular to) 
the direction of the anisotropy. 



from the boundary, the string projected onto the boundary will not be a straight line. The 
direction of the force acting on the string endpoint at the boundary can be defined by an 
angle 

tan0 = ^£. (4.13) 

If <f) = then the force acts along the z-axis. One could now think of the following 
experiment. We act with forces pointing in a specified direction in the xz-plane on two 
heavy quarks that are initially close together. We choose the strength of the forces such 
that the heavy quarks slowly start to separate more until they dissociate. When we keep the 
direction of the forces fixed the whole time the quarks will however not follow a straight 
line along the force due to the deformation of the space in the holographic coordinate. 
Instead we observe the behavior shown in Fig. 3. We note that for the JW background 
the strings bend differently depending on the sign of the B parameter. In the right panel 
of Fig. 3 we consider a string endpoint with a force acting in (j) = tt/4 direction and vary 
the depth of the turning point of the hanging string. Therefore we can probe the geometry 
up to a certain value of the holographic coordinate. As u increases the deformation of the 
string gets stronger and stronger. For strings hanging almost down to the singularity we 
notice that the strings get deformed in such a way that they smoothly fit in the remaining 
space. When we take a look at the line element of the singular gravity dual we note that 
for B > the z-direction disappears while for B < the transverse directions vanish and 
the space degenerates into an infinite line as we go to the singularity. Therefore in the JW 
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Figure 3: String profiles in the JW model for boundary forces acting on the quarks pointing in a 
fixed direction (as indicated by 4>) and projected onto the boundary. The left panel shows profiles 
for strings at the point of string breaking [E re s- = 0) for B = (black, dotted), B = y2 (blue) and 
B = —y/6 (red) with e = const. The right panel shows strings with <j) = 7r/4 but different turning 
points for the hanging string. For B = y/2 (B = —V&) the colors going from red to yellow (blue 
to green) correspond to u — |w s b, u s b, |w s b and 2u s & with it s & the value at string breaking. Here 
dashed lines indicate unstable string configurations. 
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Figure 4: Potential for heavy quarks in the MT model. The left panel compares varying 
anisotropies at constant temperature: a/T = (isotropic; black line), a/T m 1.32 (oblate; blue 
lines), a/T rj 6.43 (prolate; green) and a/T s» 50.51 (prolate; red); the right panel at constant 
entropy density for different anisotropies aNc /s 1 / 3 ~ 1.13 (oblate; blue), aN^ 3 / s 1 / 3 sw 4.23 (pro- 
late; green) and aNc^/s 1 / 3 ~ 27.37 (prolate; red). Full (dashed) lines correspond to a separation 
of the quarks along (perpendicular to) the direction of the anisotropy. 



model the pressure anisotropy is encoded very directly in the geometry which is probed by 
the hanging string. 



- 13 - 



The anisotropic plasma of the MT model dual to axion-dilaton gravity is actually in 
thermal equilibrium and therefore we can compare the heavy quark static potential at 
constant temperature and at constant entropy density. The difference can be clearly seen 
in Fig. 4. At fixed temperature the dissociation length gets smaller for any separation in 
the xz-plane as we increase the anisotropy parameter a. At constant entropy density the 
difference of the dissociation length compared to an isotropic plasma depends on whether 
we separate the quarks along a transverse direction (string breaking occurs at a larger 
distance) or along the longitudinal direction (string breaking happens at a smaller distance 
compared to the isotropic result). 

However, regardless of the sign of the pressure anisotropy A we find that in the MT 
model the heavy quark potential is always deeper for transverse separation of the quarks. 
(Note that the blue lines in the figures correspond to oblate configurations, while the green 
and red lines are for an increasingly prolate plasma, see Tables 2 and 3). This is a striking 
difference to the situation in the JW model where oblate and prolate anisotropies lead to 
opposite deformations of the heavy quark potential. The situation in the MT model is 
instead always similar to that in the JW model for prolate anisotropy. This appears to be 
a direct consequence of the fact that in the MT model g zz / 'g xx = T~L > 1 for any a whereas 
in the singular geometry of the JT model g zz j ' g xx is larger (smaller) than unity for prolate 
(oblate) pressure anisotropy. 

In the remaining plots we will only show the results for constant entropy density 
keeping in mind that for constant temperature the distance at which the string breaks 
becomes smaller and smaller as we increase a. 

Finally we also present the results for strings where the forces acting on the endpoints 
point in certain directions specified by the angle (j). In Fig. 5 we also note that the situation 
is qualitatively the same as in the prolate case for the singular gravity dual. Also indi- 
cated in the plot are the trajectories the endpoint of the string follows as we increase the 
anisotropy and keep (f> fixed. Here again we see once more that the geometry does change 
monotonically with increasing a irrespectively of the behavior of the pressure anisotropy 
in the boundary theory. In the right panel of Fig. 5 we also probe the geometry by varying 
the location of the turning point of the hanging string. 

4.2 Comparison with weak-coupling calculations 

At weak coupling, the real part of the heavy quark potential V{r) is given by the Fourier 
transform of the electrostatic propagator. In an axisymmetric situation integration over 
the azimuth angle leads to 



V{f) = -— dk dC,J Q {kr^l - C 2 sin0 r ) cos{kr( cos 6» r )L> o(^ = 0,fc), (4.14) 



where cos# r = z/r and ( = k z /k with our choice of the anisotropy direction along z. In the 
case of the hard anisotropic loops, the propagator L>oo given by (3.5) involves poles at real 
k corresponding to electric plasma instabilities which are integrated over with a principal 
value prescription, while Dqq in the zero-coupling version of the MT model only has poles 




at imaginary k. 
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Figure 5: String profiles in the MT model for boundary forces acting on the quarks pointing in 
the same direction (as indicated by </>) and projected onto the boundary. The left panel shows 
the profiles for strings at the point of string breaking [E reg - — 0) for aN^ 3 / s 1 / 3 = (black, 
dotted), aiV c 2/3 /s 1/3 ~ 1.13 (blue), aNc /3 /s^ 3 m 4.23 (green) and a7V c 2/3 /s 1/3 ~ 27.37 (red) with 
N~ 2 s « 0.159. The right panel shows the profiles for aA^/s 1 / 3 « 4.23 and 4> = tt/4 but at 
different turning points for the hanging string: at string breaking u s b (green), 4u s f,/5 (cyan) and 
6u s b/5 (purple). Here dashed lines indicate unstable string configurations. 

In Fig. 6 we have evaluated (4.14) with the hard anisotropic loop propagator for 
strongly oblate (£ = 100) and prolate (£ = —0.9) anisotropy (cf. Table 4), keeping alterna- 
tively the hard particle density n and the energy density e fixed for different anisotropies. 
The details of the deviation from the isotropic result slightly depend on whether n or e 
is kept constant, and in either case we find that for oblate anisotropies the heavy quark 
potential is slightly deeper along the anisotropy direction than transverse to it, while for 
prolate anisotropies this situation is reversed 7 . In order to make these effects more visible, 
we also plot V(L) divided by the modulus of the vacuum (Coulomb) potential, l/(47rL). 

Comparing with the results of the JW model, we find a remarkable qualitative agree- 
ment in the dependence on the sign of the anisotropy and the direction of the quark 
separation. Moreover, the absolute deviation from the isotropic result is rather small both 
at weak coupling and in the JW model. 

On the other hand, as we have seen above, the MT model has a qualitatively different 
dependence on the direction of the quark separation in the case of oblate anisotropies 
(which are usually considered in the context of heavy ion collisions). 

Turning to the zero-coupling version of the MT model introduced in Sect. 3.2, the 
heavy quark potential is given by the Fourier transform of (3.11) which is plotted in Fig. 
7. As we have discussed in Sect. 3.2, the high-temperature limit of this weak coupling 

7 This is also true when N rather than phard is rescaled in (3.1), but this method leads to somewhat 
stronger differences between fixed n and fixed e. 
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Figure 6: Static potential for a weakly coupled anisotropic plasma in the hard anisotropic loop 
formalism as a function of the quark separation L (units set by the isotropic Debye mass). The blue 
and red lines correspond to longitudinal (full) and transverse (dashed) orientation for £ = 100 and 
£ = —0.9, respectively. The isotropic result is shown in black. In the upper two plots, the particle 
number density for the anisotropic and the isotropic plasma are the same, in the lower two plots 
the energy density is kept fixed. 

model corresponds to a prolate anisotropy (in contrast to the holographic MT model), 
whereas for general a/T both prolate and oblate anisotropies are possible, depending on 
the renormalization scale. Curiously enough, the potential shown in Fig. 7 (which does 
not depend on UV renormalization) has qualitatively similar dependence on the direction 
of quark separation as the holographic MT model (and the hard anisotropic loop potential 
in the prolate case). 

We finally also consider the behavior of the quark potentials at large distances. In the 
two holographic models, there is a finite separation beyond which the string connecting 
the heavy quarks becomes unstable because strings entering the horizon or the naked 
singularity are energetically favored, and at a somewhat larger distance even no unstable 
connecting solution can be found. 

To leading order at weak coupling, the isotropic quark potential is simply given by a 
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Figure 7: Static potential in free anisotropic 0-QCD for a = 1. Full lines correspond to the 
z-direction and dashed lines to the ^-direction. 

Yukawa potential with exponential decay at large distance. The anisotropic weak coupling 
results show curious deviations. In the anisotropic ^-deformed zero-coupling case Fig. 7 
shows a nonmonotonic behavior of the potential along the anisotropy direction such that 
beyond L ~ 5 (where the potential is actually already extremely small) there is even a 
repulsive behavior. 

Even more curious behavior can be found in the hard anisotropic loop potential at 
large distances. In this case there is nonmonotonic behavior along (transverse to) the 
anisotropy direction for prolate (oblate) anisotropy, and here the nonmonotonic behavior 
is moreover oscillatory. This behavior, which has not been noted in the previous studies 
of the hard anisotropic loop potential [17, 18], is shown in Fig. 8, where the potential is 
plotted at large distances in the xz plane (enhanced by dividing by the modulus of the 
vacuum (Coulomb) potential). This oscillatory behavior, which is reminiscent of Friedel 
oscillations at finite chemical potential (for a recent discussion see [47]), has its origin in 
the presence of poles in the electrostatic propagator at real wave vector corresponding to 
electric plasma instabilities. It is however rather clear that this curious behavior is devoid 
of physical implications even at weak coupling, because the plasma instabilities imply that 
a stationary anisotropy is only a justifiable approximation at sufficiently small time scales 
and correspondingly small length scales. 

5. Jet Quenching 

5.1 Holographic calculations 

The computation of the jet quenching parameter q for an anisotropic plasma with an axion- 
dilaton-gravity dual, the MT model, has been presented in [28, 30]. Here we will reproduce 
the result for the most general case with an ultrarelativistic quark moving in an arbitrary 
direction [30] and compare the results with those of the singular geometry of the JW model. 

According to the prescription of [48, 49] we calculate the string worldsheet with end- 
points moving in the same direction at the speed of light and separated a small distance I 
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Figure 8: Ratio of static potential to vacuum potential for £ = 100 and £ = —0.9 at same energy 
density as in the isotropic plasma. 



along a direction perpendicular to their motion. The jet quenching parameter q can then 
be obtained from 



(W A (C UghtUke )) = exp (- ^=q) + (5.1) 



In the following we consider a quark endpoint moving in the xz-plane. The direction 
is parametrized by an angle 9 such that for 9 = the quark moves along the z-axis. We 
therefore start by two subsequent coordinate transformations. First we define 

Z = zcos9 + xs'mO, (5-2) 
X = x cos 9 — z sin 8, (5-3) 
Y = y (5.4) 

and then we introduce light-cone coordinates 



Z ± = -=(t±Z). (5.5) 



V2' 



The metric then takes the form 



ds 2 =G ++ {dZ + ) 2 + G—{dZ-) 2 + 2G+-dZ+dZ- (5.6) 
+ GxxdX 2 + 2G x +dXdZ + + 2G X -dXdZ~ + GyydY 2 + G vu dU 2 . 
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Writing the new metric coefficients in terms of our original ones we find 8 



C++ = G — = - (g tt + 9xx sin 2 9 + g zz cos 2 9) , (5.7) 

Gxx = 9xx cos 2 9 + g zz sin 2 9, (5.8) 

Gyy = 9xx, (5.9) 

Guu = 9uu, (5.10) 

G x + = -G X - = cos 6 sin 9(g xx - g zz ). (5.11) 



We choose the worldsheet coordinates (r, a) = (Z~ , U) and let Z + , X and Y depend 
on the holographic coordinate U in the following. It is interesting that we must allow for 
a non-constant embedding of the string in Z + to find a solution in most general case. The 
Nambu-Goto action of the string is then given by 



S = - 



1 



27TQ' 



dZ / du 



G\_{Z + ) 12 + G 2 X _X' 2 + 2G+-G X -(Z + )'X' 



G — {Guu + G ++ (Z + )' 2 + GxxX' 2 + G YY Y' 2 + 2G + x{Z + )'X') 



(5.12) 



The expression under the square root is actually negative which leads to an imaginary 
action. The reason is that we consider a spacelike string worldsheet. However this is 
expected because it is exactly what we need to obtain a jet quenching parameter that is 
real. 

Since the Lagrangian does not depend on Z + , X or Y explicitly we can find three 
constants of motion II + , Ii x and U y . In the limit where these constants are small 9 we 
obtain 



(z+) / = c++ n+ + c+x n x + o(n 2 ), (5.13) 

X' = c x +U + + cxx'nx + 0(U 2 ), (5.14) 
Y< = c Y y^y + C(n 2 ). (5.15) 



In [30] the coefficients c are given explicitly for the metric of the axion-dilaton-gravity dual. 
Since we are interested in comparing the results of two different gravity duals we express 



Here and in the following we will not write the dependence of the metric coefficients on the holographic 
coordinate u in order to keep the expressions shorter. There is no danger of confusing at which value of u we 
should evaluate the metric coefficient, because the string worldsheet in the calculation of the jet quenching 
parameter will always have its turning point at the horizon [30]. 

9 We want to consider small separation lengths I between the two string endpoints. In [30] it is shown 
that this corresponds to the limit of small II's. As a further remark we note that the worldsheet turning 
point characterized by dU/dX — (and similarly for Z + and y) is located at the horizon. 
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these coefficients in terms of the general form of the metric given in (4.1). 



9uu 9tt{gxx cos + g zz sm 0) + 

9xx9zz fr 1 c\ 

1 2(g tt + g zz cos 2 + g xx sin 2 0) 9tt9xx9z 



jzz 



2g g zz cos 2 6 + g xx sin 2 ■ 

9tt + 9zz cos 2 + # xx sin 2 g xx g z 



°xx=\ ^-z (5-17) 



CYY=\ ^ 2- — (5-18) 

V Stt + <fe cos 2 (9 + # xx sin # xx 

/ g^u (g zz - g X x) sin 0cos0 

c+x =c X + = J ; 2fl~i ( 5A9 > 

V 9tt + gzz cos^ + g xx sm z 6* fefez 

This agrees with [30] if we insert the precise form of the metric in the axion-dilaton-gravity 
case. However it is now also straightforward to consider any background whose metric is 
of the form (4.1). 

The string endpoints at the boundary are not separated along the Z + direction and 
integrating (5.13) gives 



f Q Uh du c +x 
Io h du C++ 



n+ = ~ J ° 3 + U X . (5.20) 



Along the X-axis the separation of the endpoints is I sin <j) while in Indirection it is / cos 
The constants of motion are then 

I sin 4> Jq H du C++ 

2 Q du c XX Jo h du C++ - ( fi h du c +x ) 



U x = ^ Jo "^ ++ j, (5.21) 



IV = ■ (5-22) 

2 J du cyy 

If we insert the expressions (5.13)-(5.15) into the action (5.12) and expand to second order 
in n's we obtain 

S = — duy/G—Guu + ^ du [c++n 2 + c xx U 2 x + 2c +x U + U x + c YY n Y ] . 

(5.23) 

The action is imaginary because we considered a spacelike worldsheet and L~ is the 
length of the Wilson line in ^--direction. The first, II independent term is divergent, 
however, the jet quenching parameter is proportional to I 2 and therefore just contained in 
the second finite term. Upon inserting (5.20)-(5.22) into the action and considering the 
defining relation for the jet quenching parameter (5.1) we eventually obtain 

= (P(0) sin 2 <f> + Q(0) cos 2 (j>) (5.24) 

with 

C du 

J Uh du cxx J Uh du C++ - ( J Q Uh du c +x ) 



m =— ^ ±± T^~. - 2 , (5-25) 



Q(0) =ju^ • (5-26) 

J h du C Y Y 
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Figure 9: Jet quenching parameter normalized to the isotropic result in the JW model. Left panel 
for an oblate plasma with B = v2, right panel for a prolate plasma with B = — y6, with e = const. 




Figure 10: Jet quenching parameter normalized to the isotropic result in the MT model. Left 
panel for a plasma at the same entropy density for aNc^ 3 / s 1 / 3 1.13 and N~ 2 s « 0.159 
which corresponds to oblate pressure anisotropy; right panel for prolate pressure anisotropy with 
aNc^/s 1 / 3 27.37. 

The average 

1 f 2n 1 
^ = ^ 7T J o d( P9e,<t> = ^(qe,o + qe,n/2) = Q8,n/4 (5-27) 

is the total jet quenching parameter for a quark moving with angle 9 with respect to the 
anisotropy direction, while qg^ contains the information about momentum broadening in 
directions transverse to the motion of the quark, with = being perpendicular to both 
the direction of the jet and the anisotropy direction. When 9 = 0, i.e. the quark moving 
along the direction of anisotropy, qo^ = % is independent of <fi. In the context of heavy-ion 
collisions, one is of course mostly interested in jets at larger 9. For 9 = tt/2 one can define 
transverse and longitudinal jet quenching parameters 

Q± = Qn/2,0 j QL = <?7r/2,7r/2 • (5.28) 

In Figures 9 and 10 qg^ is plotted for oblate and prolate pressure anisotropies in the 
MT model and in the JW model. Once more we see that in the JW model a different sign of 
the anisotropy parameter leads to a qualitative change of the result, in particular whether 
qi is larger or smaller than q±, while the results for the MT model are more similar to the 
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prolate phase of the JW model, for any parameter a and thus regardless of the sign of the 
pressure anisotropy A. 



5.2 Comparison with weak-coupling calculations 

Attempts to calculate the anisotropic jet quenching parameters for a quark moving trans- 
verse to the anisotropy direction {9 = vr/2) have been presented in Refs. [15, 16] based on 
one-loop calculations using hard anisotropic loops. In this situation the presence of space- 
like poles in the static gluon propagator leads to nonintegrable singularities. In Ref. [15] it 
was conjectured that these singularities might get be cured by the generation of an imagi- 
nary part in the static gluon self energy at higher loop order and this conjecture was used 
for an estimate of this "anomalous" contribution. While in Ref. [50] this conjecture was 
refuted, Ref. [16] proposed alternative resolutions, which all point to anomalous contribu- 
tions of the same sign and angular dependence as the infrared-safe regular contributions 
computed previously in Ref. [15] to leading logarithmic order. For a quark moving with 
the speed of light and in the limit of small anisotropy parameter £, the regular contribution 
to q± and qi was obtained as 



To linear order in £, the anomalous contribution is in fact only present for oblate anisotropy 
(£ > 0) with 



where C is a positive constant which depending on the physical cutoff for the singularities 
arising from plasma instabilities may differ from unity and also involve ln(l/£) [16]. At 
any rate, to linear order in £ the hard anisotropic loop calculations of Refs. [15, 16] imply 
Ql > Q± for oblate pressure anisotropy, and qi < q± for the prolate case. This result 
neither agrees with the results of the JW model nor with those of the MT model: in the 
JW model the ordering of the two jet quenching parameter changes with the sign of the 
anisotropy, but the ordering is just the opposite. The MT model on the other hand always 
has q~L > q± which agrees with the hard anisotropic loop result in the oblate case, but 
differs in the prolate case. 

It is actually questionable whether the one- loop calculation using hard anisotropic loops 
is relevant for the physics of a weakly coupled anisotropic quark-gluon plasma. Anisotropic 
jet quenching could be instead dominated by large chromomagnetic fields generated by 
plasma instabilities [51, 52]. Because plasma instabilities give rise to \B±\ > \E±\ and 
\El\ > \Bl\, it has been argued in Ref. [51] that this would also give qi > q± for a 
plasma with oblate anisotropy. In fact, the (different) plasma instabilities for prolate 
anisotropics equally lead to 11 \B±\ > \E±\ and \El\ > \Bl\, thus also favoring qi > q±. 
Perhaps fortuitously, this is in line with the results of the MT model, although the latter 

10 Refs. [15, 16] only discussed oblate anisotropies, which have m 2 a < and rr? a oc £. For prolate anisotropy 
only the last term of Eq. (32) in [16] contributes, which is of order £ 3 . 

11 This holds true for both weak fields [8] and nonperturbatively large fields [53]. 




(5.29) 




(5.30) 
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of course neglects any dynamics from instabilities and the formation of inhomogeneous 
configurations 12 . 

6. Conclusion 

In this paper we have studied two different holographic models of a strongly coupled super- 
Yang-Mills plasma with (temporarily) fixed pressure anisotropy. To this end, the JW model 
uses a singular anisotropic geometry with a naked singularity, whereas the MT model 
achieves a completely regular geometry and thus thermal equilibrium through an axion 
field with constant spacelike gradient which is dual to a parity-violating deformation of 
the gauge theory by a spatially varying 9 angle. We have probed the different geometries 
by Wilson loops which in their respective forms are dual to the static potential of heavy 
quarks and jet quenching parameters. 

In the case of the heavy quark potential, we have found that the results of the JW 
model agree qualitatively with a weak coupling description through anisotropic hard loops 
in that quarks separated moderately along the anisotropy direction have a slightly deeper 
(shallower) potential for oblate (prolate) pressure anisotropy. The MT model differs from 
the JW model in that the anisotropic deformation of the heavy quark potential is always 
similar to the prolate case of the JW model or the hard anisotropic loops, regardless of 
the sign of the pressure anisotropy. A similar result was found for a ^-deformed weakly 
coupled gauge theory which we have introduced as an alternative model for an anisotropic 
(and parity-violating) plasma. 

At larger distances in the two holographic models the quarks break up, but the disso- 
ciation lengths are ordered differently with respect to orientation and sign of the pressure 
anisotropy in oblate plasmas. On the weak coupling side we have noted a different (but 
physically rather irrelevant) behavior at large quark separation: a monotonic potential in 
some directions and oscillatory tails reminiscent of Friedel oscillations in others. 

By means of lightlike Wilson loops we have extracted anisotropic jet quenching param- 
eters for the two holographic models and compared with perturbative leading-log calcula- 
tions in hard anisotropic loops. Once again we observed qualitative changes, in particular 
of the ordering of qi and q± for transverse jets, when the pressure anisotropy is changed 
from oblate to prolate in the JW model, but not in the MT model, which is again always 
similar to the prolate case of the JW model, i.e. qi > q±. Whereas in the weak coupling 
calculations based on hard anisotropic loops the ordering of <jl and q± depends on the sign 
of the pressure anisotropy, we have argued that nonperturbatively large (inhomogeneous) 
fields generated by plasma instabilities of a weakly coupled plasma always favor q^ > q±. 

The fact that the MT model arrives at qi > q± seems to be related to the fact that 
the bulk geometry is uniformly prolate (H > 1), even when the pressure of the dual 
gauge theory has oblate anisotropy. It would be very interesting to see whether dynamical 
holographic models of heavy-ion collisions, such as colliding shock waves in AdS [22], where 
both the bulk geometry and the pressure anisotropy have oblate character, would lead to 

12 Note that the fields associated with plasma instabilities at weak coupling have nonvanishing wave 
number. 
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different results (as the results obtained in the JW model could suggest). This seems to be 
a most pertinent question since heavy-ion data [54] indeed point towards <jl > q±, which 
weak coupling calculations in an anisotropic plasma appear able to reproduce [15, 16, 51]. 

Acknowledgments 

We are indebted to Antti Gynther, Stefan Strieker, and Aleksi Vuorinen for many useful 
discussions, and to Mike Strickland for sharing computer code for the evaluation of heavy 
quark potentials in hard anisotropic loop theory. 

This work was supported by the Austrian Science Foundation FWF, project no. P22114. 

References 

[1] B. Miiller and J. L. Nagle, Results from the relativistic heavy ion collider, 
Ann.Rev.Nucl.Part.Sci. 56 (2006) 93-135, [nucl-th/0602029]. 

[2] S. Mrowczyfiski, Stream instabilities of the quark-gluon plasma, Phys. Lett. B214 (1988) 587. 

[3] P. Romatschke and M. Strickland, Collective modes of an anisotropic quark gluon plasma, 
Phys. Rev. D68 (2003) 036004, [hep-ph/0304092]. 

[4] P. Arnold, J. Lcnaghan, and G. D. Moore, QCD plasma instabilities and bottom-up 
thermalization, JEEP 08 (2003) 002, [hep-ph/0307325]. 

[5] P. B. Arnold, G. D. Moore, and L. G. Yaffc, The Fate of non-Abelian plasma instabilities in 
3+1 dimensions, Phys. Rev. D72 (2005) 054003, [hep-ph/0505212]. 

[6] A. Rebhan, P. Romatschke, and M. Strickland, Dynamics of quark- gluon-plasma instabilities 
in discretized hard-loop approximation, JEEP 0509 (2005) 041, [hep-ph/0505261]. 

[7] A. Rebhan, M. Strickland, and M. Attcms, Instabilities of an anisotropically expanding 
non-Abelian plasma: 1D+3V discretized hard-loop simulations, Phys. Rev. D78 (2008) 
045023, [0802.1714]. 

[8] A. Rebhan and D. Stcineder, Collective modes and instabilities in anisotropically expanding 
ultarelativistic plasmas, Phys. Rev. D81 (2010) 085044, [0912.5383]. 

[9] M. Attems, A. Rebhan, and M. Strickland, Instabilities of an anisotropically expanding 
non-Abelian plasma: 3D +3V discretized hard-loop simulations, 1207.5795. 

[10] P. B. Arnold and G. D. Moore, QCD plasma instabilities: The NonAbelian cascade, 
Phys. Rev. D73 (2006) 025006, [hep-ph/0509206]. 

[11] J. Berges, S. Schemer, and D. Sexty, Turbulence in nonabelian gauge theory, Phys. Lett. 
B681 (2009) 362-366, [0811.4293]. 

[12] A. Ipp, A. Rebhan, and M. Strickland, Non-Abelian plasma instabilities: SU(3) vs. SU(2), 
1012.0298. 

[13] S. Mrowczyfiski, A. Rebhan, and M. Strickland, Eard loop effective action for anisotropic 
plasmas, Phys.Rev. D70 (2004) 025004, [hep-ph/0403256]. 

[14] B. Schenke and M. Strickland, Photon production from an anisotropic quark-gluon plasma, 
Phys.Rev. D76 (2007) 025023, [hep-ph/0611332]. 



- 24 - 



P. Romatschke, Momentum broadening in an anisotropic plasma, Phys. Rev. C75 (2007) 
014901, [hep-ph/0607327]. 

R. Baier and Y. Mehtar-Tani, Jet quenching and broadening: The Transport coefficient q-hat 
in an anisotropic plasma, Phys. Rev. C78 (2008) 064906, [0806.0954]. 

A. Dumitru, Y. Guo, and M. Strickland, The Heavy-quark potential in an anisotropic 
(viscous) plasma, Phys. Lett. B662 (2008) 37-42, [0711.4722]. 

O. Philipsen and M. Tassler, On Quarkonium in an anisotropic quark gluon plasma, 
0908.1746. 

Y. Burnier, M. Laine, and M. Vepsalainen, Quarkonium dissociation in the presence of a 
small momentum space anisotropy, Phys. Lett. B678 (2009) 86-89, [0903.3467]. 

0. Aharony, S. S. Gubser, J. M. Maldacena, H. Ooguri, and Y. Oz, Large N field theories, 
string theory and gravity, Phys.Rept. 323 (2000) 183-386, [hep-th/9905111]. 

J. Casalderrey-Solana, H. Liu, D. Mateos, K. Rajagopal, and U. A. Wiedemann, 
Gauge/String Duality, Hot QCD and Heavy Ion Collisions, 1101.0618. 

P. M. Chesler and L. G. Yaffc, Holography and colliding gravitational shock waves in 
asymptotically AdS 5 spacetime, Phys. Rev. Lett. 106 (2011) 021601, [1011.3562]. 

R. A. Janik and P. Witaszczyk, Towards the description of anisotropic plasma at strong 
coupling, JHEP 0809 (2008) 026, [0806.2141]. 

A. Rebhan and D. Steineder, Electromagnetic signatures of a strongly coupled anisotropic 
plasma, JHEP 1108 (2011) 153, [1106.3539]. 

D. Mateos and D. Trancanelli, The anisotropic N=4 super Yang-Mills plasma and its 
instabilities, Phys. Rev. Lett. 107 (2011) 101601, [1105.3472]. 

D. Mateos and D. Trancanelli, Thermodynamics and Instabilities of a Strongly Coupled 
Anisotropic Plasma, JHEP 1107 (2011) 054, [1106.1637]. 

A. Rebhan and D. Steineder, Violation of the Holographic Viscosity Bound in a Strongly 
Coupled Anisotropic Plasma, Phys. Rev. Lett. 108 (2012) 021601, [1110.6825]. 

D. Giataganas, Probing strongly coupled anisotropic plasma, 1202.4436. 

M. Chcrnicoff, D. Fernandez, D. Mateos, and D. Trancanelli, Drag force in a strongly coupled 
anisotropic plasma, 1202.3696. 

M. Chernicoff, D. Fernandez, D. Mateos, and D. Trancanelli, Jet quenching in a strongly 
coupled anisotropic plasma, 1203.0561. 

K. B. Fadafan and H. Soltanpanahi, Energy loss in a strongly coupled anisotropic plasma, 
1206.2271. 

J. Erdmenger, P. Kcrner, and H. Zeller, Transport in Anisotropic Superfluids: A Holographic 
Description, JHEP 1201 (2012) 059, [1110.0007]. 

1. Gahramanov, T. Kalaydzhyan, and I. Kirsch, Anisotropic hydrodynamics, holography and 
the chiral magnetic effect, 1203.4259. 

M. Amnion, V. G. Filev, J. Tarrio, and D. Zoakos, D3/D7 Quark-Gluon Plasma with 
Magnetically Induced Anisotropy, 1207 . 1047. 



- 25 - 



[35] R. Baier, S. A. Strieker, O. Taanila, and A. Vuorincn, Holographic dilepton production in a 
thermalizing plasma, 1205.2998. 

[36] T. Azeyanagi, W. Li, and T. Takayanagi, On String Theory Duals of Lifshitz-like Fixed 
Points, JEEP 0906 (2009) 084, [0905.0688]. 

[37] E. Braaten and R. D. Pisarski, Simple effective Lagrangian for hard thermal loops, Phys. Rev. 
D45 (1992) 1827-1830. 

[38] M. Margotta, K. McCarty, C. McGahan, M. Strickland, and D. Yager-Elorriaga, Quarkonium 
states in a complex-valued potential, Phys. Rev. D83 (2011) 105019, [1101.4651]. 

[39] W.-T. Ni, Equivalence Principles and Electromagnetism, Phys. Rev. Lett. 38 (1977) 301-304. 

[40] J. Ralston, Gauge theory in the adiabatic approximation, Phys. Rev. D51 (1995) 2018-2021. 

[41] S. M. Carroll, G. B. Field, and R. Jackiw, Limits on a Lorentz and Parity Violating 
Modification of Electrodynamics, Phys. Rev. D41 (1990) 1231. 

[42] A. Gynther, A. Rebhan, and D. Steineder, Thermodynamics and phase diagram of 
anisotropic Chern-Simons deformed gauge theories, 1207 . 6283. 

[43] J. M. Maldacena, Wilson loops in large N field theories, Phys. Rev. Lett. 80 (1998) 4859-4862, 
[hep-th/9803002]. 

[44] S.-J. Rey, S. Theisen, and J.-T. Yee, Wilson-Polyakov loop at finite temperature in large N 
gauge theory and anti-de Sitter supergravity, Nucl.Phys. B527 (1998) 171-186, 
[hep-th/9803135]. 

[45] A. Brandhuber, N. Itzhaki, J. Sonnenschein, and S. Yankiclowicz, Wilson loops in the large N 
limit at finite temperature, Phys. Lett. B434 (1998) 36-40, [hep-th/9803137]. 

[46] M. Laine, O. Philipscn, P. Romatschke, and M. Tassler, Real-time static potential in hot 
QCD, JEEP 0703 (2007) 054, [hep-ph/0611300]. 

[47] J. -P. Blaizot, A. Ipp, and A. Rebhan, Study of the gluon propagator in the large-Nf limit at 
finite temperature and chemical potential for weak and strong couplings, Annals Phys. 321 
(2006) 2128-2155, [hep-ph/0508317]. 

[48] H. Liu, K. Rajagopal, and U. A. Wiedemann, Calculating the jet quenching parameter from 
AdS/CFT, Phys. Rev. Lett. 97 (2006) 182301, [riep-ph/0605178]. 

[49] H. Liu, K. Rajagopal, and U. A. Wiedemann, Wilson loops in heavy ion collisions and their 
calculation in AdS/CFT, JEEP 0703 (2007) 066, [hep-ph/0612168]. 

[50] M. Carrington and A. Rebhan, On the imaginary part of the next-to-leading- order static 
gluon self-energy in an anisotropic plasma, Phys. Rev. D80 (2009) 065035, [0906.5220]. 

[51] A. Dumitru, Y. Nara, B. Schenke, and M. Strickland, Jet broadening in unstable non-Abelian 
plasmas, Phys. Rev. C78 (2008) 024909, [0710.1223]. 

[52] A. Majumder, B. Miiller, and S. Mrowczyhski, Momentum Broadening of a Fast Parton in a 
Perturbative Quark-Gluon Plasma, Phys. Rev. D80 (2009) 125020, [0903.3683]. 

[53] A. Ipp. private communication. 

[54] J. Putschke, Lntra-jet correlations of high-p t hadrons from STAR, J. Phys. G 34 (2007) 
S679-684, [nucl-ex/0701074]. 



- 26 - 



